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We study the effects on the spectrum and distribution of high-energy neutrinos due to scattering 
with dark matter both outside and within our galaxy, focusing on the neutrinos observed by the 
IceCube experiment with energies up to several PeV. If these neutrinos originate from extra-galactic 
astrophysical sources, then scattering in transit with dark matter particles will delay their arrival 
to Earth. This results in a cut-off in their spectrum at an energy set by the scattering cross section, 
allowing us to place an upper limit on cross sections cr which increase with energy E at the level of 
a < • (m/GeV) • (E/PeV)^ cm^, for dark matter particles of mass m. Once these neutrinos 

enter our galaxy, the large dark matter densities result in further scattering, especially towards the 
Galactic Centre. Intriguingly, we find that for cr ~ • (m/GeV)• (E/PeV)^ cm^, the distribution 

of the neutrinos on the sky has a small cluster of events towards the centre of the galaxy, potentially 
explaining the ~ 2 sigma excess seen by IceCube in this region without needing a galactic source. 


I. INTRODUCTION 

Dark matter (DM) is inferred to exist, through its grav¬ 
itational interactions with visible matter, within and be¬ 
tween galaxies [1-4]. However its nature at the parti¬ 
cle level will remain unknown unless its interactions with 
the particles of the Standard Model are understood. One 
of the least constrained interactions of DM is with neu¬ 
trinos. Bounds exist on the cross section for scattering 
at low energy from the Cosmic Microwave Background 
(CMB) and large-scale structure [5-7] and limits also ex¬ 
ist at MeV energies from supernovae [6-9]. However the 
DM-v scattering cross section is only weakly constrained 
at higher neutrino energies (TeV-scale and larger) due to 
a relative lack of observational tests, both astrophysical 
and terrestrial. Such high-energy probes are important 
as they provide information on how DM particles interact 
at energies inaccessible to colliders, for example. 

We show that the recent observation of PeV-energy 
neutrinos by the IceCube experiment [10] allows the 
DM-v scattering cross section to be probed at much 
higher energies than has previously been possible (see 
also refs. [11-15] for other new physics probes). In¬ 
deed the IceCube experiment has now observed 37 neu¬ 
trino events with energies between 30 and 2000 TeV after 
988 days of data-taking. The distribution of these neutri¬ 
nos is approximately isotropic, implying an extra-galactic 
origin (see also ref. [16]), however there is potentially a 
small cluster of events towards the Galactic Centre. 

In order to study the interactions of these high-energy 
neutrinos with DM, we consider two different regimes: 
scattering of neutrinos with DM between galaxies and 
within the Milky Way itself. For the former we show in 
section 11 that a large scattering cross section results in a 


delayed arrival to our galaxy of the neutrinos, since they 
must diffuse towards Earth instead of free-streaming. 
Hence the observation of PeV-energy neutrinos by Ice- 
Cube means that they have not been impeded signifi¬ 
cantly in transit. By making the assumption that the 
highest energy events observed by IceCube originate from 
extra-galactic hypernova remnants [17-19], we calculate 
the cross section a for DM-z/ scattering which would re¬ 
sult in a cut-off in the observed spectrum below 2 PeV, 
allowing us to set an upper limit on a. 

For the case of scattering within our galaxy we show 
in section HI that DM-z/ scattering results in the genera¬ 
tion of an anisotropic distribution of neutrino events from 
a purely extra-galactic flux, due to enhanced scattering 
within the Galactic Centre region where the DM density 
is largest. We calculate the cross section needed to gen¬ 
erate the amount of clustering of high-energy neutrino 
events towards the Galactic Centre potentially observed 
by the IceCube experiment [10]. Finally in section IV we 
summarise our findings and conclude. 


II. DARK MATTER-NEUTRINO SCATTERING 
IN THE INTERGALACTIC MEDIUM 

Due to their weak interactions with ordinary mat¬ 
ter, the PeV neutrinos observed by IceCube should have 
travelled unimpeded from their extra-galactic sources to 
Earth. However this may not be the case as the rate 
at which these neutrinos scatter with dark matter par¬ 
ticles is poorly constrained [5-8] and could be large, re¬ 
sulting in them diffusing towards our galaxy instead of 
free-streaming. In this section, we constrain the scatter¬ 
ing cross section a for DM between galaxies using the 
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FIG. 1: Transit of a PeV-energy neutrino from its extra- 
galactic source to Earth at two redshift values 5. With weak 
DM-i/ scattering, the path of the neutrino is lengthened be¬ 
tween the source and Earth, causing it to arrive at the present 
day 5 = 0. When the scattering rate is large, i.e. strong scat¬ 
tering, the neutrino has to diffuse towards Earth, in which 
case it has not arrived by the present day. 


spectrum of high-energy neutrinos observed by IceCube. 

The mean free path for scattering is Adm-i^ = 
w/(pdmO'), where m is the dark matter mass and pdm 
is its density in the intergalactic medium (IGM). We as¬ 
sume that the latter takes the form pdm = /0o(l + 
where po = 2.7 • 10“^^ kgm“^ and z is the redshift [20]. 
If Adm-i/ is smaller than the distance between galaxies, 
then the neutrinos will scatter multiple times, length¬ 
ening their path from the source to Earth and thereby 
delaying their arrival. We assume that the DM mass 
is much larger than that of the neutrinos and that the 
scattering is elastic and flavour-universal [21]. 

A schematic of this effect is shown in figure 1, where 
we give three broad regimes of scattering. In the case of 
‘weak scattering’, we have that a > 0 and so the path the 
neutrino has to travel is elongated, causing it to arrive 
later at the present day (z = 0). The case of ‘strong 
scattering’ where cr ^ 0 has the neutrino scattering so 
often as it travels from outside our galaxy that it has yet 
to arrive at Earth by the present day. 

We now use this formalism to calculate the effect on 
the flux of extra-galactic neutrinos. Without scattering 
between DM and neutrinos, the expected flux is obtained 
by integrating over sources at all redshifts z [17-19, 22], 



dz R{z) dA^ 

47r77o v/iiM(lTzp~TDA 


Ki + z)E], ( 1 ) 


where dN/dE is the neutrino spectrum at the source, 
R{z) is the rate of the particular source events for the 
neutrinos, Hq = 0.69 kms“^Mpc“^ is the Hubble con¬ 
stant, Dm = 0.27 is the fraction of matter in the Universe 
and Da = 0.73 is the fraction of dark energy [20]. We 
assume that R{z) is proportional to the star formation 
rate i?sFR(-2:) [17-19]. 

Since the effect of DM-neutrino scattering is to delay 


the arrival of the neutrinos to Earth, the redshift at which 
the neutrino is emitted is no longer necessarily equal to 
the redshift of the source. Hence the integral over redshift 
from equation (1) is split into two integrals: one over the 
redshift distance to the source z' and another over the 
emission time of the neutrinos z, 



R{z)f(.z',z,a) dN 
•\/Dm(1 + z)^ + Da ‘^E 


[{l + z)E]- 


( 2 ) 

The time delay due to the DM-neutrino scattering is 
contained in the (normalised) delay function /(z',i, cr), 
which gives the probability that a neutrino was emitted 
at a redshift z, given that its source is at a distance corre¬ 
sponding to redshift z'. The amount of delay depends on 
the cross section of scattering between dark matter and 
neutrinos cr, which could itself vary with the neutrino en¬ 
ergy E. For figure 2 and throughout this work we assume 
that a = ao{m/GeV) ■ (E/Eq) , where Eq = 1 PeV. This 
is an empirical parameterisation, and to preserve unitar- 
ity, the cross section should decrease at higher energies 
above a few PeV. Hence in this case f{z',z,a) also de¬ 
pends on the neutrino energy as it travels to Earth. 

Examples of f{z', z, a) which have been calculated nu¬ 
merically for different values of uq are shown in fig¬ 
ure 2, assuming a source distance of 10 Mpc correspond¬ 
ing to z' « 0.0025 and a neutrino energy at the source 
of I PeV. For erg = 10“^® cm^ ,the neutrinos reach¬ 
ing Earth have been produced recently, and so values 
of z which are close to z' will contribute to the inte¬ 
gral in equation (2). Indeed when (Tq = 0 iWe have that 
f{z',z,a = 0) = 5{z' — z), since the neutrino was emit¬ 
ted at the same redshift as we observe the source to be 
at today and eq. (2) reduces to eq. (1). 

By contrast for erg = 10“^® cm^ and trg = 5-10“^® cm^, 
the longer path means the neutrinos originated from their 
source farther in the past. Hence for these cross sections, 
the integral in equation (2) is dominated by values of z 



FIG. 2: Fraction of neutrinos f{z', z, a) which were produced 
at a given time in the past represented by a value of redshift 
z, given that they originate from an extra-galactic source at 
a distance of 10 Mpc from Earth (at the present day). Larger 
values of the DM-i^ scattering cross section a mean that neu¬ 
trinos arriving today were produced longer ago. We also show 
the star formation rate as a function of z for comparison. 
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FIG. 3: Spectra of high-energy neutrinos from extra-galactic 
hypernova remnants for different values of their scattering 
cross section with DM. These are compared with the case of 
zero scattering and data from the IceCube experiment [10]. 
We assume a cross section which increases with neutrino en¬ 
ergy E as a — cro(m/GeV) • (i?/PeV)^. Larger cross sections 
result in a cut-off in the spectrum at lower energies. 


much larger than z'. As such for ctq > 5 • 10“^® cm^, 
the majority of neutrinos arriving at the present day will 
have been produced at 5 > 5, when the star-formation 
rate was small. In this case, the neutrino flux would likely 
be much smaller than that observed by IceCube. 

We now perform the full integral in equation 2 by 
assuming a specific astrophysical model for the high- 
energy neutrino flux observed by IceCube. We use the 
model proposed in refs. [17-19] in which the highest- 
energy (> 100 TeV) neutrinos originate from the de¬ 
cay of hadrons accelerated at extra-galactic hypernova 
remnants (HNR) (alternative models for these neutri¬ 
nos use e.g. Active Galactic Nuclei [23] or gamma- 
ray bursts [24]). We calculate the expected flux using 
eqn. (2), replacing R{z) with the HNR rate which we 
take to be R{z) = 10 ~‘^Mq^Rsfr{z) [17-19]. For the 
source neutrino spectrum dN/dE, we use the formalism 
of refs. [17-19], based on a convolution of the spectrum of 
protons accelerated at the HNR with the secondary neu¬ 
trino spectrum produced by each proton (via e.g. pion 
or muon decay). This depends on the maximum energy 
to which these protons can be accelerated, which we take 
to be = IQi® eV for the HNR [17-19]. 

The effect of scattering between DM and high-energy 
neutrinos in transit to Earth from their extra-galactic 
HNR sources is shown in figure 3, compared with the as¬ 
trophysical neutrino spectrum observed by IceCube [10, 
25]. The result is a cut-off in the spectrum of neutrinos 
at an energy set by the size of (Tq. This arises from neu¬ 
trinos which scatter so often that they do not arrive at 
Earth by the present day, even if they are produced when 
the first hypernovae form around 5 ~ 5. Larger values 
of tTo effectively widen the range of energies for which 
this condition holds. Since as can be seen in figure 3 
the IceCube data extend up to A ~ 2 PeV we conserva¬ 
tively set an upper limit on a DM-iz cross section which 
increases with at the level of (Tq ^ 10“^^ cm^. This 


limit is robust to changes in the model for the production 
of these neutrinos, as the cut-off is a generic feature of 
DM-i/ scattering. Future measurements of higher energy 
neutrinos will serve to strengthen this upper limit. 

III. SCATTERING WITHIN THE MILKY WAY 

In this section, we consider the effects of high-energy 
extra-galactic neutrinos scattering with DM within our 
own galaxy, after they have travelled through the IGM. 
We have already shown that DM-i/ scattering between 
galaxies results in a cut-off in the spectrum of these neu¬ 
trinos for large enough cross sections. However here we 
show that scattering within the Milky Way, particularly 
at the centre of the galaxy, affects the distribution of 
high-energy neutrino events on the sky. This is interest¬ 
ing as the IceCube collaboration have reported a small 
(;< 2(7 significance) clustering of high-energy neutrino 
events towards the Galactic Centre [10, 26], which corre¬ 
sponds to approximately 15% of all of their events. 

We assume that the neutrinos enter the galaxy from all 
directions isotropically as expected if their origin is extra- 
galactic. However even a kiloparsec from the Galac¬ 
tic Centre, the dark matter density can reach values as 
large as p ~ 100 GeVcm“^ for a Navarro-Frenk-White 
(NFW) profile [27], and the density increases further to¬ 
wards the centre. Hence at this DM density, a value 
of (To > 10“^^ cm^ for a PeV-energy neutrino would be 
enough for the mean free path Adm-i/ to be of kiloparsec 
scale or smaller, and so there can be significant scattering 
between DM and neutrinos towards the Galactic Centre 
even given the constraints of the previous section. This 
will lead to neutrinos appearing to originate from this 
region even if they are produced outside of our galaxy. 

The effect of DM-z/ scattering within our galaxy on the 
distribution of extra-galactic neutrinos (with an energy 
of 1 PeV) is shown in figure 4 for two different values 
of the scattering cross section (assuming an NFW profile 
for the DM). For (Tq = 10“^^ cm^, the scattering between 
DM and PeV-energy neutrinos generates a concentration 
of events towards the Galactic Centre, where the DM 
density (and therefore the scattering rate) is largest, at 
the level of about 15 — 20% of the total rate. In this case 
scattering has altered the direction of travel for these neu¬ 
trinos and therefore their apparent point of origin. By 
contrast for (Tq = 10“^^ cm^ (and below), the neutrino 
flux remains mostly isotropic. Hence (Tq = 10“^^ cm^ 
gives a similar amount of clustering towards the Galactic 
Centre seen in the latest IceCube data [10, 26], while val¬ 
ues larger than this would introduce too much anisotropy. 

The Galactic Centre is not the only region for which 
an anisotropy could be observed: the large DM densities 
towards the cores of dwarf galaxies [28-30] could result 
in a similar effect from DM-z^ scattering. The size of 
this anisotropy depends on the mass model for the dwarf 
galaxies, which is not well known. For models where 
the DM density is > 100 GeVcm“^ within 1 kpc of the 
centre [29], the mean free path with ao = 10“^^ cm^ is 
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FIG. 4: Distribution of neutrinos (with E = 1 PeV) on the sky centred on the Galactic Centre for two different values of 
the scattering cross section with DM cr = cro • (m/GeV). The (logarithmic) colour scale gives the fraction of the total rate in 
each angular pixel. We assume that the neutrinos are produced outside the galaxy and enter isotropically, however for DM-v 
scattering within our galaxy with ag = 10“^^ cm^ the high densities of DM in the Galactic Centre region create a cluster of 
events. For ag = cm^ by contrast the scattering rate is much lower and so the neutrino flux remains mostly isotropic. 


smaller than 1 kpc. Hence the flux of high energy neutri¬ 
nos around e.g. the Fornax or Sculptor dwarfs could be 
enhanced by as much as a few percent compared to the 
average all-sky flux, depending on the density profile. 

IV. CONCLUSION 

We have investigated the effects of scattering between 
high-energy neutrinos and dark matter particles, focusing 
on the neutrino events observed by the IceCube experi¬ 
ment between energies of £1 « 30 — 2000 TeV [10]. Since 
there are strong arguments for an extra-galactic origin 
of these neutrinos [10, 16] we have studied the effect of 
scattering both outside and within our galaxy. For extra- 
galactic DM-z/ scattering with a cross section which rises 
with we showed that this results in a delayed arrival 
for the neutrinos to Earth (figures 1 and 2), due to the 
large distances (> 10 Mpc) they must diffuse through 
when the scattering rate is large. This leads to a cut¬ 
off in their spectrum at an energy set by the size of the 
cross section (figure 3). While for scattering within our 
galaxy, we found that the distribution of the neutrinos 
is shifted from pure isotropy to possessing an excess of 
events towards the Galactic Centre, where the DM den¬ 
sity is largest (figure 4). This results from scattering al¬ 
tering the direction of travel for these neutrinos, making 
them appear to originate from the Galactic Centre. 



a{m/GeV) ^ [cm^] 

Energy-scale 

This work (gal.) 
This work (ex-gal) 
Structure form. [5] 
SN 1987A [8] 

~ 10-52 . (E/eV)^ 

< 10-4^ • {E/eY f 

< 10-58 • {E/eWf 

< 10-28 

100 - 2000 TeV 

100 - 2000 TeV 

< eV 

~ 10 MeV 


TABLE I: Here ‘gal.’ refers to best-fit cross section cr to 
give the required level of clustering seen towards the Galactic 
Centre in the recent IceCube data [10]. ‘Ex-gal’ is our upper 
limit from requiring that neutrinos with energies of 2 PeV and 
below are able to reach Earth by the present day. 


We summarise our numerical results compared with 
those from previous works in table I. As can be seen, our 
results probe much higher energy scales than previous 
constraints from structure formation [5, 6] and supernova 
1987A [6-8]. Throughout this work, we have assumed 
that the scattering cross section a increases as up to 
a few PeV (or above), since in this case our constraint is 
particularly strong compared to those set at lower ener¬ 
gies. However this is not the only form which the DM-ic 
scattering cross section can take, and the complimentar- 
ily between the various constraints serves as an important 
probe of the energy dependence of cr. The upper bound 
labelled ‘ex-gal’ in table I has been set by assuming that 
the neutrinos originate from hypernova remnants in other 
galaxies, although any other sources for these neutrinos, 
astrophysical [23, 24] or otherwise [26, 31-34], would lead 
to similar constraints. For cross sections larger than this 
upper limit, neutrinos with energies below 2 PeV are im¬ 
peded too much on their way to our galaxy, which would 
lead to a cut-off in the IceCube neutrino spectrum at this 
energy. Scattering cross sections just below this bound 
give no observable cut-off while also leading to interest¬ 
ing phenomenology within our galaxy. Indeed the best-fit 
cross section to give a small cluster of PeV-energy events 
(^ 15 — 20%) towards the Galactic Centre, due to DM-i^ 
scattering within our galaxy, is ^ 10“^^ • (m/GeV) cm^ 
(labelled ‘gal’ in table I). This is consistent with the < 2cr 
anisotropy observed in the IceCube data-set in this re¬ 
gion [10, 26]. Cross-sections of this size which scale with 
E'^ also evade the low-energy constraints in table I. 

Currently there are not enough statistics to make a 
more precise statement. However we predict that if DM 
particles scatter with ~PeV energy neutrinos with a cross 
section around 10“^^ • (m/CeV) cm^, then future data 
from IceCube (and experiments such as KM3NeT [35] 
or Antares [36]) should show neutrinos with a mostly 
isotropic distribution plus an excess of events (~ 15% of 
the population) towards the Galactic Centre. The neutri¬ 
nos are produced in extra-galactic astrophysical sources, 
with the DM only contributing to the anisotropy through 
scattering, and so the spectrum will retain the general 
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power-law forni expected from such sources regardless of 
position on the sky (although the exact spectrum de¬ 
pends on the energy dependence of the cross section). 
This combined spectral and spatial dependence is a po¬ 
tentially unique prediction compared to alternative mod¬ 
els for the anisotropy. For example although an astro- 
physical source within our galaxy (on top of an extra- 
galactic component) may result in a similar prediction, it 
would likely over-produce other signals (e.g. 7 or cosmic 
rays) [16, 37]. Additionally results from Antares dis¬ 
favour a point source [36]. Furthermore models where the 
neutrinos are produced by decaying or annihilating DM 
lead to an excess of neutrinos towards the Galactic Cen¬ 
tre, but the spectrum of these events will be very different 
from those away from this region [26, 31-34]. We also ex¬ 
pect similar ^ 10 % anisotropies towards other galaxies 


such as M87, within the regions they subtend on the sky, 
and potentially small anisotropies for the dwarf galaxies 
near the Milky Way. Our results are complimentary to 
bounds on DM-^ interactions at lower energies [5-8] and 
also to probes of scattering of e.g. DM with quarks [38- 
40], DM with electrons [41], DM with photons [42, 43] 
and self-interactions [44], placing us closer than ever to 
understanding the nature of DM. 
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